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Abstract. We compute Zero Point Energy in a spherically symmetric background with the help 
of the Wheeler-DeWitt equation. This last one is regarded as a Sturm-Liouville problem with the 
cosmological constant considered as the associated eigenvalue. The graviton contribution, at one 
loop is extracted with the help of a variational approach together with Gaussian trial functionals. 
The divergences handled with a zeta function regularization are compared with the results obtained 
using a Noncommutative Geometry (NCG) and Modified Dispersion Relations (MDR). In both 
NCG and MDR no renormalization scheme is necessary to remove infinities in contrast to what 
happens in conventional approaches. Effects on photon propagation are briefly discussed. 
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INTRODUCTION 

The Cosmological Constant problem is certainly one of the most fascinating challenges 
of our days. A challenge because all the attempts that try to explain the 10^^^ orders of 
magnitude of discrepancy between the theory and observation have produced unsatisfy- 
ing results. If we believe that Quantum Field Theory (QFT) is a part of the real world, 
then the theoretical predictive power to compute the Cosmological Constant must be 
entrusted to the methods of QFT at the Planck scale. Indeed, calculating the Zero Point 
Energy (ZPE) of some field of mass m with a cutoff at the Planck scale, we obtain 

EzPE = \ v/^^ ^ ^ \Q^'GeV\ (1) 

while the observation leads to a ZPE of the order 10^^^ GeV^ . This surely represents 
one of the worst predictions of QFT. However if one insists to use the methods of 
QFT applied to General Relativity, one necessarily meets one of the most famous 
equations appeared in the literature of Cosmology and Gravity: the Wheeler-DeWitt 
(WDW) equation[l]. The WDW equation was originally introduced by Bryce DeWitt as 
an attempt to quantize General Relativity in a Hamiltonian formulation. It is described 



by 



(2K)Gijki7t'j7t''-^{'R-2A) 



(2) 



and it represents the quantum version of the classical constraint which guarantees the 
invariance under time reparametrization. Gijki is the super-metric, Tt'-' is the super- 
momentum,^i? is the scalar curvature in three dimensions and A is the cosmological 
constant, while K = SnG with G the Newton's constant. An immediate application of 
the WDW Equation is given in terms of the Friedmann-Robertson-Walker (FRW) mini 
superspace, where all the degrees of freedom but the scale factor are frozen. The FRW 
metric is described by the following line element 



ds'^ = -N^dt^ + a'^ {t)dal, 



where dQ.^ is the usual line element on the three sphere, namely 



d^l = Yijdx^dxK 



In this background, we have simply 
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and the WDW equation Hm (a) = 0, becomes 



da da 
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Eq.(6) assumes the familiar form of a one-dimensional Schrodinger equation for a 
particle moving in the potential 



U{a) = 
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with total zero energy. The parameter q represents the factor-ordering ambiguity and 

ao = is a reference length. The WDW equation (6) has been solved exactly in 
terms of Airy functions by Vilenkin[2] for the special case of operator ordering q = 
— 1. The Cosmological Constant A here appears as a parameter. Nevertheless, except 
the FRW case and other few examples, the WDW equation is very difficult to solve. 
This difficulty increases considerably when the mini superspace approach is avoided. 
However some information can be gained if one changes the point of view. Indeed, 
instead of treating A in Eq.(2) as a parameter, one can formally rewrite the WDW 
equation as an expectation value computation^ [3]. Indeed, if we multiply Eq.(2) by 



See also Ref.[4] for an application of the method to a / (R) theory. 



^* [gij] and functionally integrate over the three spatial metric gij we find 



1 J^[gij] [gtj] f^d^xA^^ [gtj] _ 1 (^|/^j3^Az|^) _ A 

In Eq.(8) we have also integrated over the hypersurface L and we have defined 

V^j^d\^ (9) 

as the volume of the hypersurface L with 

As = {2k) Gijkin'^n^^ - ^^R/ {2k) . (10) 

In this form, Eq.(8) can be used to compute ZPE provided that A/jc be considered as 

an eigenvalue of A^. In particular, Eq.(8) represents the Sturm-Liouville problem as- 
sociated with the cosmological constant. To solve Eq.(8) is a quite impossible task. 
Therefore, we are oriented to use a variational approach with trial wave functional. 
The related boundary conditions are dictated by the choice of the trial wave function- 
als which, in our case are of the Gaussian type. Different types of wave functionals 
correspond to different boundary conditions. The choice of a Gaussian wave functional 
is justified by the fact that ZPE should be described by a good candidate of the ''vac- 
uum state". In the next section we give the general guidelines in ordinary gravity and in 
presence of Modified Dispersion Relations and the Non Commutative approach to QFT. 
Units in which ^ = c = fc=lare used throughout the paper. 



HIGH ENERGY GRAVITY MODIFICATION: THE EXAMPLE OF 
NON COMMUTATIVE THEORIES AND GRAVITY'S RAINBOW 



As an application of the Eq. (8), we consider the simple example of the Mini Superspace 
described by a FRW cosmology (3). We find the following simple expectation value 



JS'a'¥* {a) 



4G^ 



/^a^* {a) [fl4]^(a) 



3A;r^ 



(11) 



where the normalization is modified by a weight factor. The application of a variational 
procedure with a trial wave functional of the form 



^ = exp {-pa^) 



(12) 



shows that there is no real solution of the parameter j8 compatible with the procedure. 
Nevertheless, a couple of imaginary solutions of the variational parameter 



^ 4G 



(13) 



can be found. These solutions could be related to the tunneling wave function of 
Vilenkin[2]. Even if the eigenvalue procedure of Eq.(ll) leads to an imaginary cos- 
mological constant, which does not correspond to the measured observable, it does not 
mean that the procedure is useless. Indeed, Eq.(8) can be used to calculate the cosmo- 
logical constant induced by quantum fluctuations of the gravitational field. To fix ideas, 
we choose the following form of the metric 



ds^ = -N^{r)dt'^ + 



1 

r 



+ ? {dd'^ + singed 



(14) 



where b (r) is subject to the only condition b{rt) — r^. We consider gij = gij + where 
gij is the background metric and hij is a quantum fluctuation around the background. 
Then we expand Eq.(8) in terms of hij. Since the kinetic part of Aj: is quadratic in the 
momenta, we only need to expand the three- scalar curvature / d^x^^R up to second 
order in hij^. As shown in Ref.[6], the final contribution does not include ghosts and 
simply becomes 
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The integration over Gaussian Wave functionals leads to 



where 



Arh' 



(16) 



(17) 



is the modified Lichnerowicz operator and A/,is the Lichnerowicz operator defined by 



(18) 



Qijki represents the inverse DeWitt metric and all indices run from one to three. Note that 
the term —AR\h^--[- ^Rhjj disappears in four dimensions. The propagator (^;^)/aW 
be represented as 

where h^^^^ (^) are the eigenfunctions of A^. T denotes a complete set of indices and 
A (t) are a set of variational parameters to be determined by the minimization of Eq. ( 1 6) . 
The expectation value of A^ is easily obtained by inserting the form of the propagator 



^ See Refs.[4, 5, 6] for technical details. 



into Eq.( 16) and minimizing with respect to the variational function A (t). Thus the total 
one loop energy density for TT tensors becomes 



A 
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The above expression makes sense only for cof (t) > 0, where Ct); are the eigenvalues of 
A^. Following Refs.[4, 5, 6], we find that the final evaluation of expression (20) is 
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where we have included an additional An coming from the angular integration and where 
we have defined two r-dependent effective masses m\ (r) and m\ (r) 



^l{r)^j,(\-'^)+^b'{r) + ^b{r) 



{r = r{x)). 



(22) 



The effective masses have different expression from case to case. For example, in the 
Schwarzschild case, b (r) = 2MG, we find 
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(23) 



The expression in Eq.(21) is divergent and must be regularized. For example, the zeta 
function regularization method leads to 



e 



mf (r) ^/eJ _ 



/= 1,2 



(24) 



where an additional mass parameter ji has been introduced in order to restore the 
correct dimension for the regularized quantities. Such an arbitrary mass scale emerges 
unavoidably in any regularization scheme. The renormalization is performed via the 
absorption of the divergent part into the re-definition of a bare classical quantity and the 
final result is given by 



Ao _ ml (r) 
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Of course, one can follow other methods to obtain finite results: for instance, the use 
of a UV-cut off. Nevertheless it is possible to obtain finite results introducing a distor- 
tion in the space-time from the beginning. This can be realized with the help of the 
Non Commutative Approach to QFT developed in Ref.[7] or with the help of Gravity's 



Rainbow developed in Ref.[9]. Noncommutative theories provide a powerful method to 
naturally regularize divergent integrals appearing in QFT. Eq.(21) is a typical example of 
a divergent integral. The noncommutativity of spacetime is encoded in the commutator 
[x'^, x^] — i 6^^, where 0^"^ is an antisymmetric matrix which determines the fundamen- 
tal discretization of spacetime. In even dimensional space-time, 6^^ can be brought to a 
block-diagonal form by a suitable Lorentz rotation leading to 



e^" = diag 01 e"^ 028"^ . . . 0rf/2e 



^ab 



(26) 



with e"^ a 2 X 2 antisymmetric Ricci Levi-Civita tensor. If 0, = Vz = 1 . . .(i/2, then the 
space-time is homogeneous and preserves isotropy. The effect of the length on ZPE 
calculation is basically the following: the classical Liouville counting number of nodes 



dn 



(27) 



is modified by distorting the counting of nodes in the following way [7] 

d^xd^k . d^xd^k 
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This deformation corresponds to an effective cut off on the background geometry (14). 
The UV cut off is triggered only by higher momenta modes > 1/ \/0 which propagate 
over the background geometry. As an effect the final induced cosmological constant 
becomes 
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where an integration by parts in Eq.(21) has been done. We recover the usual divergent 
integral when — )^ 0. The result is finite and we have an induced cosmological constant 
which is regular. We can obtain enough information in the asymptotic regimes when the 
background satisfies the relation 



ml (r) = m\ (r) = -ml (r) , 



(30) 



which is valid for the Schwarzschild, Schwarzschild-de Sitter (SdS) and Schwarzschild- 
Anti de Sitter (SAdS) metric close to the throat. Indeed, defining 



X = 
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(31) 



we find that when x -|-oo, 
A 



^TtG en 
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Conversely, when x ^ 0, we obtain 
A 4 
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The other interesting cases, namely de Sitter and Anti-de Sitter and Minkowski are 
described by 



mj (r) = m2 (r) = mo (r) , 



leading to 
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when X — )■ 0. As regards Gravity's Rainbow[8], we can begin by defining a'^rainbow 
metric" 



+ ■ 



{dd^ + sm^ed<j)^). (37) 



8iiE/Ep) ' (^i_bM.^g2^E/Ep) 8l{E/Ep) 
gi (E/Ep) and g2 (E/Ep) are two arbitrary functions which have the following property 



lim gi{E/Ep)^l 



and lim g2(E/Ep) — l. 

E/Ep-^0 
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We expect the functions gi {E/Ep) and g2 {E/Ep) modify the UV behavior in the same 
way as GUP and Noncommutative geometry do, respectively. Following Ref.[9], in 
presence of Gravity's Rainbow, we find that Eq.(8) changes into 
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Of course, Eq. (39) and Eq. (40) reduce to the ordinary Eqs. (2, 8) and (10) when E /Ep — > 
0. By repeating the procedure leading to Eq. (20), we find that the TT tensor contribution 
of Eq.(39) to the total one loop energy density becomes 
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where E* is the value which annihilates the argument of the root. In the previous 
equation we have assumed that the effective mass does not depend on the energy E. 
To further proceed, we choose a form of gi (E/Ep) and g2{E/Ep) which allows a 
comparison with the results obtained with a Noncommutative geometry computation 
expressed by Eq.(29). We are thus led to choose 

gi{E/Ep) = ^M-OC^) and giiE/Ep) = 1, (42) 

with a > and j8 e M, because the pure ''Gaussian" choice with j8 = can not give a 
positive induced cosmological constant^. However this is true when the effective masses 
satisfy relation (30). In case relation (34) holds the pure "Gaussian" choice works 

for large and small x, where x = itiq (r) /Ep. The final result is a vanishing induced 
cosmological constant in both asymptotic regimes. It is interesting to note that Gravity's 
Rainbow has potential effects on the photon propagation[10]. Indeed, let us consider two 
photons emitted at the same time t — —to at x^s — 0- The first photon be a low energy 
photon (E <^ Ep) and the second one be a Planckian photon (E ~ Ep). Both photons are 
assumed to be detected at a later time in x^s- We expect to detect the two photons with a 
time delay At given by the solution of the equation 

xf<<^-(0)=4-^-(A^), (43) 

that implies 

\/\ff/^ 

A? ^ gl {E) ^/X^ ^^Tp^'V^ \l^eff/^to) , (44) 

where we have used (42) for the rainbow functions. 
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